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The couphng to gravity in D=5 spacetime dimensions is considered for the 
particle-hke and vortex-type solutions obtained by uplifting the D=4 Yang- 
Mills instantons and D=3 Yang-Mills-Higgs monopoles. It turns out that the 
particles become completely destroyed by gravity, while the vortices admit 
a rich spectrum of gravitating generalizations. Such vortex defects may be 
interesting in view of the AdS/CFT correspondence or in the context of the 
brane world scenario. 
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1. Introduction.— The recent interest in the AdS/CFT correspondence 
(see [|l| for a review), and in the brane world scenario |10| has attracted 
attention to gravity theories in D=5 spacetime dimensions. In this connec- 
tion solutions of gauged supergravities have been actively studied. As such 
theories generically contain the non-Abelian gauge fields, considering solu- 
tions for such fields coupled to gravity seems to be important. At the same 
time, most studies so far have been restricted exclusively to the Abelian sec- 
tor. Practically all what is known about gravitating non-Abelian solutions 
in D=5 are the BPS configurations described in In D=4 on the other 
hand, gravitating Yang-Mills (YM) fields have been analyzed in some detail 
(see for a review), partially with the numerical methods. 

The aim of this letter is to study more systematically gravity-coupled 
YM fields in D=5 by applying the methods previously used in D=4. Instead 
of specializing to any particular supergravity model, we shall consider the 
pure Einstein- YM (EYM) theory. Although this theory is probably non- 
super symmetric in itself, it enters all gauged supergravities as the basic 
building block, and hence one can expect some features of its solutions to be 
generic. In principle, one can study this theory also in higher dimensions, 
but we shall work only in D=5. 

As is well known, the pure YM theory in D=5 Minkowski space admits 
topologically stable, particle-like and vortex-type solutions obtained by up- 
lifting the D=4 YM instantons and D=3 YM-Higgs monopoles. It is then 
natural to wonder what happens to these objects when they are coupled to 
gravity. This issue seems to be interesting in itself, and it has not been ad- 
dressed before. Below we shall study this problem and find that the particles 
become completely destroyed by gravity, as a result of their scaling behavior, 
while the vortices admit very non-trivial gravitating generalizations. These 
gravitating vortices comprise an infinite family including the fundamental 
solution and its excitations. All solutions exist for a finite range of values 
of the gravitational coupling constant, and in the strong gravity limit they 
become gravitationally closed. 

It seems plausible that such solutions could be further generalized to 
include scalars and a cosmological term. In this case they would describe 
stable vortex excitations over the AdSs background. Such topological defects 
might be interesting in view of the AdS/CFT correspondence as providing 
the dual gravity description for some processes on the gauge theory side, or 
possibly in the context of the brane world scenario. 
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2. The model.— The pure Einstein- YM theory for the gauge group SU(2) 
is defined by the action 



Here F^,^ = OmA^ - d^Al, + SabcAljAj, (a = 1, 2, 3), and [G'/^] = [g^] = 



[length]. Let us split the coordinates as = {x^,x^), where = 
In the flat spacetime limit, G 0, the YM theory is not conformally in- 
variant, the length scale being (7^, and this allows for soliton solutions. 
Specifically, for static, purely magnetic fields with = (0,A^(x'^)) the 
energy E = ^ J(F^j^yd'^x coincides with the action of the D=4 Euclidean 
YM theory. It follows then that there are regular, topologically stable solu- 
tions in D=5 with the energy E = 8'n''^\n\/ g"^, where the topological wind- 
ing number n G 7r3(SU(2)). These solutions describe neutral, particle-like 
objects which we shall call "YM instanton particles". If d/dx'^ is a sym- 
metry, one can choose A%i = (0, and the energy per unit 
x^ E = ^Jiid^H" + eabcA\H^f + \{F^,;)'^)d^x, coincides with the energy 
of the D=3 YM-Higgs fields. Since H'^H'^ is asymptotically constant [§], 
there are topological vortex-type solutions with the energy per unit length 
E = 4:n\n\/g'^, where n G 712(5"^). These we shall call "YM vortices". Let us 
now set G 7^ 0. 

3. Gravitating YM particles.— We parameterize the S0(4)-invariant line 
element in the one-particle sector as 



with = 1 — Km/r^, the length scale being g"^. The gauge field is given 
by A" = (1 + w) 6°". Here a, m are functions of r, rffig = O'^O'^ is the 
line element of the unit sphere 5*^, and 9*^ are the invariant forms on 5*^, 
such that dO"- + SabcO'^ A 6''^ = 0. The gravitational coupling constant is 
K = 87rG/g^. The ADM mass, M = m(oo), determines the dimensionful 
energy -Eadm = (3vr^/(yf^)M. The independent field equations read 




(1) 



ds' = a'Ndt^- — -r'dnl 



(2) 



r'^Nw" + rw' + K(m 
rm' = r'^Nw''^ + {w"^ 



r 



2 {w^ - l)w , 



(3) 
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Figure 1: Typical solution to Eqs.(||) (with k, = 10 ^) 

while the equation for a decouples and can be integrated giving cr(r) = 
exp (^—K ^ dr^. It follows that r{may = [r'^w''^ + (w^ — l)^)cr, which can 
be integrated with the boundary conditions m(0) = to give 

M[w, k]= j^'^ [rV + [w^ - If) exp (^-k ^ . (4) 

This is the ADM mass for configurations subject to the (00) and (rr) Einstein 
equations, and for such fields it is proportional to the action density. It follows 
that on-shell fields are stationary points of M[w, k]. For k = the integrand 
in @) can be represented as a sum of a total derivative and a perfect square. 
The latter vanishes if the self-duality equation holds, rw' = ±(w^ — 1), whose 
solutions are w = t 112^\ , where 6 is a scale parameter. These describe the 
regular "instanton particles" with the energy -Eadm = Svr^/gf^. 

It turns out that for k 7^ Eqs. (|^) do not admit globally regular solutions 
with finite ADM mass. Indeed, for such solutions a{r) 7^ 0, such that the 
exponent in (P never vanishes, and therefore (^ can only be finite if w{r) — >■ 
±1 for r 0, 00. It follows then that M[w{Xr), k] = M[w{r), A^k], where A 
is a constant scaling parameter. Since M should be stationary with respect 
to rescalings, one should have = -^M[w{r), \'^k]\\=i, but the latter is 
manifestly negative for k 7^ 0. 
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The YM particles therefore do not generahze to curved space - as a con- 
sequence of their scahng behavior. Specifically, since the size of the particles 
can be arbitrary, there is the corresponding scaling zero mode of the energy 
functional. This does not allow one to apply the inverse function theorem 
[^] to construct the gravitating generalizations even for arbitrarily small k. 
Qualitatively, the YM particles resemble dust, since their rescaling changes 
the gravitational binding (the exponential factor in (^), without affecting 
the energy of the gauge field (the rest of the integrand in (^)). As a result, 
the attraction and repulsion are not balanced and the equilibrium is impos- 
sible - the gravitating system wants to shrink but there is no pressure to 
stop the contraction. Notice that these arguments apply only in D=5. In 
D=4 the YM instantons can be coupled to gravity without problems - since 
their energy-momentum tensor then vanishes, they fulfill the coupled system 
of the Euclidean EYM equations for any value of G. 

One can numerically integrate Eqs.(0) with the regular boundary con- 
ditions at the origin, w = 1 — 26r^ + O(r^), m = O(r^), in order to see 
what actually happens to the static YM particle solution if k ^ 0. It turns 
out that as r increases, w first follows very closely (for small n) the flat 
space configuration, but then it fails to reach the value — 1 as r — oo and 
starts replicating the same pattern infinitely many times. As a result there 
emerges the quasi-periodic structure shown in Fig.|l]. This can be explained: 
integrating the Yang-Mills equation in (^ gives 

w'^ - {w^ - 1)' = - r (ln{a^N)yw'^dT , (5) 

with ' = = r\fN ^. This describes a particle moving with friction in the 
inverted double- well potential U = — (w^ — 1)^. At r = — oo the particle starts 
at the local maximum of the potential at w = 1, but then it looses energy 
due to the dissipation and cannot reach the second maximum at w = — 1. 
Hence it bounces back. For large r the dissipative term tends to zero, and 
the particle ends up oscillating in the potential well with a constant energy. 
After each oscillation the mass function m increases in a step-like fashion, 
and for large r one has m ~ r ~ Inr. As a result, there emerges an infinite 
sequence of static spherical shells of the YM energy in the D=5 spacetime. 
4. Gravitating YM vortices.— We parameterize the S0(3)-invariant met- 
ric in the one- vortex sector as 

ds^ = e^-dt' - e'^dr^ - d^l - e^^ {dx')\ (6) 
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with = d'd'^ + sin^ {^dip"^. The S0(3)-invariant gauge field is given by 

A^ = {w-1) Eabc n^drf + Hn'' dx^ , (7) 

where = (sint? cos sint? sin yj, cos-;?). Here z/, A, /i, w, and H are 
functions of r. With s = e'^+(+'^t^-^ the field equations read 

s (s e-^f^wj - e^'+^'^iw^ - l)w = e^^+^'^H^w , (8) 
s(se-2^/7')' = 2e2^+2A.^2^^ (g) 

s (s u')' ='^{A + 2B + 2C + D), (10) 

s(sC')'= ^(-2A + 25-4C + D), (11) 

s (s /i')' - e2'^+2^+2^ = ^ (A - - C - 2D) , (12) 

+ 2z/V' + 2Cy + C'l^') - 6^"^^^+^^ 
= ^{A + 2B -2C - D) . (13) 

Here A = e-^'^s'^H'^, B = q-^^s^w'^, C = e'^'+^^w'^H^ , D = e^^+^'^{w^ - 1)^ 
and ' = ^. In this system the last equation is the initial value "Gauss" 
constraint generating the residual gauge symmetry r f(r). This symmetry 
can be used to impose a gauge condition on the amplitudes. The equations 
are also invariant with respect to the rescalings of the t and coordinates, 
u ^ u + Uq, ( ^ ( + (q, H ^ e'^°H, with constant and Co- 
in addition, the equations are invariant with respect to the dilatations, 
fi ^ fi + e, X~^X + e, (-^C + e, k, —>■ e^^n, with constant e. The associated 
conserved Noether charge is 

Q = s{-p' + C + KHH'e-^^). (14) 

The origin of this symmetry can be traced to the fact that effectively the 
system can be viewed as the D=4 EYM-Higgs theory coupled to a dilaton. 
In general, if nothing depends on x^, one can parameterize the 5- fields as 
gMNdx^dx^ = e-^'^^.ydx^'dx'' - e^'^{dx^f and Aj^dx^ = A^dx^" + H'^dx^. 
The action (|I|) then reduces to = / dx'^ J d^x^ £4 with 

£4 = ^ + ^ ^^^^^^ - i ^^F;^F^^^ + 1 e-^'^V^H^V^H^ . (15) 
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Here TZ is the Ricci scalar for 7^1,, the indices are hfted by 7'^'^, also = 
df,Al - dyAI + SabcA^Al and V^.H" = df,H^ + eabcA^^H''. This determines the 
theory of coupled EYM-Higgs-dilaton fields in D=4. The amplitude C, there- 
fore effectively plays the role of the dilaton, which explains the dilatational 
symmetry. The vortex solutions under consideration thus can be viewed as 
D=4 gravitating YM-Higgs monopoles coupled to the dilaton in a special 
way (the systems considered so far in the literature do not have the direct 
coupling between the dilaton and Higgs fields p.) 

Some simplest solutions of the system (|8|)- ([T3D can be found. The vacuum 
Schwarzschild solution is obtained for w = ±1, H = 0, e~^^ = 1 — and 
either u = —A, ( = or ( = —A, z/ = 0. The extreme charged Abelian black 
string solution is described hj w = 0, H = Ho, e^'^ = e^*" = e""*" = 1 — rn/r, 

with rn = ^J2k/3. 

We are interested in globally regular solutions with finite mass per unit 
x^. This implies that z/ = A = C = Oatr = 00, and if e'^ = r then 
u = ~M + 0(r-2) for large r, while A = |:Madm + 0{r-^). Here M and 
^ADM, respectively, are the Newtonian and the ADM masses, the dimen- 
sionful energy being E = [Att / g'^)M . Choosing e'^ = r, dividing Eq.(|TU]) by 
{ks) and integrating gives 

M= rdre" (-e-^-^H'' + e-^w" + e^-^w''H^ + e^+< (w^ - ^^^^ 
Jo \2 2r^ J 

In fiat space, with u = X = ( = 0, the integrand here can be rearranged as 
the sum of a total derivative and two perfect squares. The latter vanish if 
the Bogomol'nyi equations hold, r^H' + w"^ — 1 = and w' + wH = 0. The 
BPS solution, H = cothr — 1/r and w = r/ sinhr, describes the fiat space 
vortex with M = Madm = 1- 

The Noether charge Q vanishes for globally regular solutions. As a result, 
using (|Hp with Q = allows us to exclude z/ from the equations. Next, 
introducing h = e~^H the amplitude ( can also be eliminated. Choosing the 
radial gauge where A = and defining R = e^^, the independent variables 
are w, h, Z = and R. If these are determined, the amplitudes u and ( 
are given by C = / Zdr, and u = ( + k JiZh"^ + hh')dr. The field equations 
(PD-(p!2D can be reformulated as a seven-dimensional autonomous system 

-^yk = Fk{ys,K), (17) 
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with yk = {w, w', h, h', Z, R, R'}, where the functions Fkiys, i^) can be read- 
off from In addition, the constraint (13) will restrict the initial 
values of the y^'s. Eqs.(|T^) possess the dilatational symmetry 

r — > er, w—^w, h^h/e, Z^Z/e, R ^ eR, K^e^K, (18) 

which also changes the mass as M — M/e. The following fixed points of the 
equations will be important in our analysis: 

I. The origin, {w, h, Z, R) = (1,0,0,0). Here the constraint "stable 
manifold" of solutions that are regular for r ^ can be characterized by the 
Taylor expansions 

w = l-6r^ + 0(r^), h = ar + 0{r^), 

Z = 0{r^), R = r + 0{r^). (19) 

II. Infinity, {w, h, Z, R) = (0,f,0,oo), where the (unconstrained) stable 
manifold can be described by 

w = Ar^^e"^^ + o(e-"), Z = nQr'^ + 0{r~^ Inr), 

h = v{l-Cr-^) + 0{r-'^\xir), (20) 

i? = r — m In r + rri^r"^ Inr — roo + "^r^^ + 0(r^^ Inr) . 

Here m = k{Cv'^ + (2 + kv'^)Q) determines the ADM mass, iV/ADM — ~ 
while M = IMadm - fQ- In (|r9D,(|20|) a, 6, A, v, C, Q, r^, 7 are eight free 
parameters, but due to the scaling symmetry ([T8| ) only seven are independent, 
and we shall set f = 1. 

HI. "Warped" AdSs x S"^. This fixed point is determined by the real root 
of the algebraic equation + + llg = 1: 

,,2_„ p2_ .2_ (1-9) (V-13g+l) 

^' (4g2-13g+l)' i?2 ' ^ {4q+l)R^ ' 

(21) 

Evaluating, 

1 27 31 

y; = 0.29, h = -^, Z = ±-^, R=0.75y^. (22) 

This is a new exact, essentially non-Abelian solution to the field equations. 
Its geometry is 

ds^ = e^^^^de - dr^ - e^^^ (dx^ - R^ dnj, (23) 
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Figure 2: The vortex solution in the strong gravity regime. 

with a = 1 + k/i^ = 2.62 (notice that for a = 1 this would be the metric 
on 74^5*3 X 5*^). The gauge field is given by with H = e^^'h. Linearizing 
Eqs.(|T7p around this solution, one finds the 7 characteristic eigenvalues of 
the linearized system to be {-2.77, -2.47, -2.12, -0.61 ± 1.24 i, 0.88, 1.54}, 
in units of l/y/R. The solution ( PT| ) is therefore a hyperbolic (saddle) fixed 
point with five stable and two unstable for r ^ oo eigenmodes. 

We now integrate the equations ( p!7| ) for a given k 7^ in the interval 
r G [0,00) using (p!^) , (pO|) as the boundary conditions and adjusting the 
seven independent free parameters to match the seven i/kS. This gives the 
gravitating vortex solutions. For small k they are only slightly deformed 
as compared to the fiat space vortex. The essentially non-Abelian field is 
contained in the central core where w 0, while in the far field zone w 
vanishes and only the long range U(l) component of the gauge field survives. 
The metric amplitudes are Zf» 0, i?~r, i?'i=al. 

It turns out to be convenient to measure the coupling to gravity not by 
K, but by ^ = OK where a comes from ([Tol) . Solutions are specified by their 
asymptotic parameters in (|T9|) , (pO|) , which lie on the curve r(^) in the eight- 
dimensional parameter space with coordinates F = {K,,a,b, A,C,Q,r^,'j). 
As ^ increases and F(^) deviates from F(0), the solutions deviate more and 
more from their fiat space profiles: R'{r) develops more and more pronounced 
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Figure 3: The exterior limiting solution. 



minimum at some finite r where Z{r) considerably deviates from zero; see 
Fig.|^. It is interesting that k not always increases with ^, but only up to 
the maximal value Kmax = 3.22 and then starts decreasing; see Fig.^ where 
projections of r(^) on the Ck, CQ, CA, and C7 planes are shown. 

For strongly gravitating solutions, as ^ continues to increase, the ampli- 
tudes w{r), h{r), Z{r), R{r) start developing oscillations around the constant 
values which are close to those given by Eq.(^). The number of oscillations 
grows as gravity gets stronger, while their amplitudes tend to zero. The 
proper length of the interval where oscillations take place increases, and so- 
lutions develop a long throat connecting the core of the vortex with the 
asymptotic region. In the strong gravity limit this throat gets infinitely long, 
and the vortex core becomes disconnected from the outside world. A similar 
phenomenon has been observed in D=4 p, 

In order to qualitatively understand such a behavior, each solution can 
be viewed as a trajectory interpolating in the phase space between the two 
fixed points (|^) and (pOl) . On its way, it gets attracted by the saddle point 
(pTD, due to the five stable for r — 00 eigenmodes around this point, it 
spends some "time" in its vicinity, but finally it gets repelled due to the two 
unstable eigenmodes. As gravity gets stronger with growing ^, the trajec- 
tory approaches closer and closer the saddle point (^I]) oscillating longer and 
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Figure 4: Parameters of the gravitating vortex solutions. The right ends 
of the curves correspond to the flat space vortex, ^ = 0. As increases, 
the curves extend to the left, and for large ^ they spiral towards the values 
corresponding to the limiting solution. 

longer in its vicinity. Finally, the trajectory exactly hits the saddle point 
and splits into two parts. The first part starts at the origin (]T^) and af- 



ter infinitely many oscillations arrives at the saddle point (|21|). The second 
part interpolates between the saddle point and infinity (^^. As a result, 
the vortex solution splits in the limit into two independent solutions. The 
first, interior solution contains the regular core, but asymptotically it is not 
flat and approaches instead the geometry (pSf). The second, exterior solution 
interpolates between (p3D and the asymptotically flat region. 

One can directly construct the exterior limiting solution. This lives in the 
interval r G {—oo, oo). Shifting r to set Too = leaves four free coefficients in 
(POD, while the (constraint) set of solutions that approach for r —>■ — cxd 
is determined by the two unstable modes around this fixed point. There are 
altogether six free parameters, and the matching conditions can be fulfilled 
if only K, is treated as the seventh free parameter. The solution is found for 
K, = 0.316; see Fig.^ This solution can be viewed as an extreme non-Abelian 
black string: in r = R coordinates e~^^ has a double zero at r/j = 0.42, while 



e ~ 



. , /'^ _ 7-^)2.02 g2C U r/j)°-'^'^ for r Vh- 
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It is instructive to consider the behavior of the parameters of the solutions 
(k, a, 6, A, C, Q, Too, 7) with increasing ^. One can shift r — r — Too to set 
'"oo = 0, which changes A Ae'^°°, 7—^7 + mroo without affecting the 
other parameters. One can then plot all parameters against ^, but it is more 
illustrative to pick up one parameter, C, say, and plot the remaining ones 
against it. The resulting curves (some of them are shown in Fig.^]) exhibit 
the characteristic spiraling behavior in the strong gravity limit (large ^). A 
similar oscillatory behavior of parameters in the critical limit is known for 
relativistic/boson stars 0. 

Apart from the fundamental solutions, one finds also their excitations, for 
which w oscillates around zero. These can be parameterized by the number 
of nodes n = l,2,...ofw. These solutions do not have the flat space limit: 
for K — > their mass M ~ 1/ This limit can be analyzed with the 
rescaling (|T^ with e = which eliminates the explicit K-dependence 

from the equations. The rescaled Higgs field then decouples in the limit, and 
this leads to solutions which resemble the D=4 Bartnik-McKinnon solutions 
[0] - with w oscillating around zero and w{oo) = dzl. 
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